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Abstract — A class of two-bit bit flipping algorithms for decod- 
ing low-density parity-check codes over the binary symmetric 
channel was proposed in 1 1 1. Initial results showed that decoders 
which employ a group of these algorithms operating in parallel 
can offer low error floor decoding for high-speed applications. As 
the number of two-bit bit flipping algorithms is large, designing 
such a decoder is not a trivial task. In this paper, we describe 
a procedure to select collections of algorithms that work well 
together. This procedure relies on a recursive process which 
enumerates error configurations that are uncorrectable by a given 
algorithm. The error configurations uncorrectable by a given 
algorithm form its trapping set profile. Based on their trapping 
set profiles, algorithms are selected so that in parallel, they can 
correct a fixed number of errors with high probability. 

I. Introduction 

With the introduction of high speed applications such as 
flash memory, fiber and free-space optical communications 
comes the need for fast and low-complexity error control 
coding. Message passing algorithms for decoding low-density 
parity-check (LDPC) codes such as the sum-product algorithm 
(SPA) offer very attractive error performance, especially for 
codes with column-weight d c > 4. However, the complexity of 
these algorithms is still high and the decoding speed is limited, 
mostly due the fact that the operations at variable and check 
nodes must be carried out for every edge in the Tanner graph. 
For regular column-weight-three LDPC codes, which allow 
lower complexity implementation, message passing algorithms 
(as well as other classes of decoding algorithms) usually suffer 
from high error floor. This weakness of message passing al- 
gorithms in regular column-weight-three LDPC codes justifies 
the search for alternatives which offer better trade-offs between 
complexity, decoding speed and error performance. 

Among existing decoding algorithms for LDPC codes on 
the binary symmetric channel (BSC), bit flipping algorithms 
are the fastest and least complex. The check node operations 
of these algorithms are modulo-two additions while the vari- 
able node operations are simple comparisons. The simplicity 
of these algorithms also makes them amenable to analysis. 
Many important and interesting results on the error correction 
capability of the serial and parallel bit flipping algorithms have 
been derived (see JTJ for a list of references). Unfortunately, 
their error performance is typically inferior. As a result, bit- 
flipping-oriented algorithms have been largely considered to 
be impractical, even after the introduction of some improved 
versions, such as the one in |(2). 



In (TJ, a class of bit flipping algorithms that employ two bits 
for decoding LDPC codes over the BSC was proposed. Com- 
pared to serial and parallel bit flipping, a two-bit bit flipping 
algorithm employs one additional bit at a variable node and 
one at a check node. The additional bits introduce memory in 
the decoding process, which slows down the decoding when 
necessary. Initial results showed that decoders which employ a 
group of these algorithms operating in parallel lower the error 
floor while maintaining low complexity. However, in |l) we 
have not given a complete failure analysis of these algorithms, 
nor have we established the methodology to derive good 
algorithms and/or a collection of mutually good algorithms. 

In this paper, we provide complete failure analysis for 
two-bit bit flipping algorithms. More importantly, we give a 
rigorous procedure to select groups of algorithms based on 
their complementariness in correcting different error patterns. 
Decoders that employ algorithms selected using this procedure 
have provably good error performance and, by the nature of bit 
flipping, high speed. As a result, two-bit bit flipping algorithms 
can finally be considered practical and suitable for high speed 
applications. 

As one can expect, a two-bit bit flipping algorithm (like 
other sub-optimal graph-decoding algorithms) fails on some 
low-weight error patterns due to the presence of certain small 
subgraphs in the Tanner graph. In this paper, we characterize 
a special class of these subgraphs and refer to them with the 
common term "trapping sets." Our definition of a trapping 
set for a given algorithm readily gives a sufficient condition 
for successful decoding. The set of all possible trapping sets 
of a given decoding algorithm constitutes the algorithm's 
trapping set profile. A unique property of trapping sets for 
two-bit bit flipping algorithms is that a trapping set profile may 
be obtained by a recursive procedure. The diversity among 
trapping set profiles of different algorithms allows us to select 
groups of algorithms such that they can collectively correct 
enor patterns that are uncorrectable by individual algorithms. 

The rest of the paper is organized as follows. Section|il]gives 
the necessary background and reviews the class of two-bit bit 



flipping algorithms. Section III gives motivation. In Section 
|IV| we define trapping sets, trapping set profiles and describe 
the recursive procedure for constructing a trapping set pro- 
file. Section [V] discusses the process of selecting algorithms. 
Numerical results are given in Section VI 



II. Preliminaries 

Let C denote an (n, k) binary LDPC code. C is defined by 
the null space of H, an m x n parity-check matrix of C. H is 
the bi-adjacency matrix of G, a Tanner graph representation 
of C. G is a bipartite graph with two sets of nodes: n 
variable (bit) nodes V(G) = {1,2,..., n} and m check nodes 
C{G) = {1,2,..., m}; and a set of edges E(G). A (d Vl d c )- 
regular LDPC code has a Tanner graph G in which all variable 
nodes have degree d v and all check nodes have degree d c . In 
this paper, we only consider (d v , d c ) -regular LDPC codes. A 
subgraph of a bipartite graph G is a bipartite graph U such 
that V(U) C V(G), C(U) C C(G) and E{U) C E(G). 
G is said to contain U, Furthermore, if Y is a graph which 
is isomorphic to U then G is also said to contain Y. In a 
bipartite graph G, the induced subgraph on a set of variable 
nodes V s C V(G'), is a bipartite graph £7 with V(U) = V s , 
C{U) = {c G C(G) : 3v G V s such that (v, c) G and 
£(*7) = {(«, c) G £(G) : u G V s }. 

A vector x = (xi,X2, . . . ,x n ) is a codeword if and only 
if xH T — 0, where H T is the transpose of H. Assume 
the transmission of the all-zero code word over the BSC. 
Denote by y the channel output vector and denote by x' = 
{x\,x l 2 i ■ ■ ■ , x n ) the decision vector after the I th iteration of 
the iterative algorithm, where I is a positive integer. At the end 
of the Z th iteration, a variable node v is said to be corrupt if 
x v — 1, otherwise it is correct. For the sake of convenience, we 
let x° = y. A variable node v with x® = 1 is initially corrupt, 
otherwise it is initially correct. Let s l — (s x , s 2 , . . . , s m ) 
denote the syndrome vector of the decision vector after the 
I th iteration, i.e., s' = x. l H T . A check node c is said to be 
satisfied at the beginning of the I th iteration if s^T 1 = 0, 
otherwise it is unsatisfied. 

Two-bit bit flipping algorithms are defined as follows. 

Definition 1: The class & of two-bit bit flipping algorithms 
is given in Algorithm[T] where the vector z l = (z[, z 2 , ■ ■ ■ , z l m ) 
gives the states of the check nodes at the beginning of 
the I th iteration while the vector w' = (w[ , w 2 , ■ ■ ■ , w n ) 
gives the states of the variable nodes at the end of the 
I th iteration. A variable node v takes its state from the set 
A v = {Os, W , l w , l s }, i.e., it can be strong zero, weak zero, 
weak one or strong one. A check node takes its state from the 
set A c = {0 p , n , l p , l n }, i.e., it can be previously satisfied, 
newly satisfied, previously unsatisfied or newly unsatisfied. 
The state wJJ of a variable node v is initialized to A v (0) £ 
{0 S ,0 W } if y v = and to A v (l) e {1 B , l w } if y v = 1. The 
state w\ of a check node c is initialized to A c (0) € {0 P , n } 
if s ° = and to A c (l) G {l p ,l n } otherwise. A two- 
bit bit flipping algorithm T = (/, Z™,A V ,A C ) iteratively 
updates z l and w' until all check nodes are satisfied or 
until a maximum number of iteration is reached. The 
check node update function $ : {0, l} 2 — > A c is defined 
as follows: $(0,0) = P ,$(0,1) = l n ,$(l,0) = n and 
$(1, 1) = l p . The variable node update is specified by a 
function / : A v x — > A v , where E.d v is the set of 
all ordered 4-tuples £, = (6,^2,^3,^4) sucn that & G N 



and Ei^ = d v- In Algorithm [T Xo p (. v )>Xo n (v), Xi p (v) 
and Xi a ( v ) gi ye tne number of check nodes with states 
z c = p ,0 n ,l p and l n , respectively, that are connected to 
v. The function / must be symmetric with respect to and 1 
and must allow every state of a variable node to be reachable 
from any other state. 

Algorithm 1 Two-Bit Bit Flipping Algorithm 

Vv: w°+- A v (y v ), Vc: zl 4- A c (s°), I «- 1 
while s l ^ and I < l™ do 

W : w l v <r- f(w l -\ X l o < (v),X l o n ( v )iXi p (v),Xi a ( v ))' 
Vc : z l c +1 <P{s l c -\s f c ); 
I 4-1 + 1; 
end while 

What makes a two-bit bit flipping algorithm novel is that a 
variable nodes has "strength" and a check node's reliability is 
evaluated based on its state in the previous iteration. 

III. Motivation 

Consider a collection stf of iterative decoding algorithms 
for LDPC codes. Let us assume for a moment that the set of 
all uncorrectable error patterns for each and every algorithms 
in is known. More precisely, in the context of LDPC 
codes, we assume that the induced subgraphs on such error 
patterns can be enumerated for each decoding algorithm. This 
naturally suggests the use of a decoder T> which employs 
multiple algorithms drawn from s/ . The basis for this use of 
multiple algorithms is rather simple: If different algorithms are 
capable of correcting different error patterns, then a decoder 
employing a set of properly selected algorithms can achieve 
provably better error performance than any single-algorithm 
decoder. Disappointingly, the above hypothetical assumption 
is not valid for most iterative algorithms. For message passing 
algorithms such as the SPA, there is no simple criterion to 
verify weather or not an arbitrary error pattern is correctable, 
much less an explicit methodology to design a decoder which 
employs multiple algorithms in a collaborative manner. 

Interestingly, for two-bit bit flipping algorithms, we are 
able to establish a framework to analyze and enumerate all 
uncorrectable error patterns, and this is the main contribution 
of this paper. In particular, we characterize the decoding fail- 
ures of two-bit bit flipping algorithms by redefining trapping 
sets and introducing the definition of trapping set profiles. 
It is an important property of the newly defined trapping 
sets that enable us to enumerate them using a recursive 
procedure. We remark that the enumeration of trapping sets is 
code independent. More importantly, the concept and explicit 
construction of trapping set profiles allow rigorous selections 
of multiple algorithms which can collectively correct a fixed 
number of errors with high probability. Given that the selection 
of multiple algorithms would become straightforward once 
the trapping sets/trapping set profiles have been defined and 
constructed, we devote a considerable portion of the paper to 
introducing these two objects. We also focus on giving criteria 



for selecting algorithms rather than explicitly describing the 
selection process. 

IV. Trapping Sets and Trapping Set Profiles 

A. Trapping Sets of Two-Bit Bit Flipping Algorithms 

Although the term trapping set was originally defined as a 
set of variable nodes that are not eventually correctable by an 
iterative decoding algorithm [3], in the literature it has been 
used more frequently to refer to a combinatorially defined 
subgraph that may be harmful to decoding. The justification 
for this less rigorous use of terminology is that the variable 
node set of a so-called trapping set (a subgraph) would be an 
actual set of non-eventually-correctable variable nodes if the 
parallel bit flipping algorithm were used (see |4| for details). 
Such a trapping set is a subgraph S in which all variable 
nodes are connected to more even-degree check nodes than 
odd-degree check nodes and no two odd-degree check nodes 
share a variable node v ^ V(S) |4J. It is important to note 
that these subgraphs are responsible for the most frequent (but 
not all) failures of the parallel bit flipping algorithm. Our 
analysis reveals that besides the above-mentioned subgraphs, 
other subgraphs in which variable nodes are connected to more 
check nodes with degree p > 1 than those with degree p = 1 
can also be harmful. Because two-bit bit flipping algorithms 
in general have a stronger error correction capability than the 
parallel bit flipping algorithm, it is necessary to include all 
harmful subgraphs into the analysis. To achieve this goal, we 
therefore (re)define the notion of a trapping set for two-bit 
bit flipping algorithms, as we now explain. We first introduce 
the following definition on failures of a two-bit bit flipping 
algorithm. 

Definition 2: Consider a two-bit bit flipping algorithm T = 
(/, A v , A c ) and a Tanner graph G. Let V c denote the set 
of variable nodes that are initially corrupt and let / denote 
the induced subgraph on V e . If the algorithm T does not 
converge on G after I™ iterations, then we say that T fails 
on the subgraph I of G. 

It can be seen that the decoding failure of T is defined 
with the knowledge of the induced subgraph on the set of 
initially corrupt variable nodes. To characterize failures of J 7 , 
a collection of all induced subgraphs / must be enumerated. 
While this is difficult in general, for practically important cases 
of small numbers of initial errors (less than 10) and small 
column-weight codes (d v = 3 or 4), the enumeration of such 
induced subgraphs is tractable. 

Consider a given Tanner graph /. Let S'i(J-) denote a set of 
Tanner graphs containing a subgraph J isomorphic to / such 
that T fails on J. Since Sj (J 7 ) is undeniably too general to be 
useful, we focus our attention on a subset S'f(J-) of <£/(.F), 
described as follows. 

Let Si <E S'liJ 7 ) such that T fails on the subgraph J\ of 
Si. Then, Si G S'i'iJ 7 ) if there does not exist S 2 € Si(T) 
such that: (i) T fails on the subgraph J2 of S%, and (ii) there 
is an isomorphism between S2 and a subgraph of Si under 
which the variable node set V(J2) is mapped bijectively into 
the variable node set V(Ji). 



Now we are ready to define trapping sets and trapping set 
profiles of a two-bit bit flipping algorithm. 

Definition 3: If S £ ^/(-T 7 ) then S is a trapping set of J- '. 
I is called an inducing set of S. S'j(J-) is called the trapping 
set profile with inducing set I of T . 

The following proposition states an important property of a 
trapping set. 

Proposition 1: Let S be a trapping set of T with inducing 
set /. Consider the decoding of T on S with V(I) being the 
set of initially corrupt variable nodes. Then, for any variable 
node v e V(S), there exist an integer < I < Z™ such that 

w l v e {i s , i„}. 

Proof: The proof is omitted due to page limits. ■ 
From Proposition [T] one can see that the trapping set profile 
Sj (J 7 ) of T contains the graphs that are most "compact." We 
consider these graphs most compact because the decoding of 
T on such a graph could be made successful by removing 
any variable node of the graph. An important remark is that 
although we do not have a combinatorial definition of trapping 
sets, there exists an explicit recursive procedure to obtain all 
trapping sets up to a certain size. We also remark that for 
certain reasonably good algorithms, the necessary condition 
for a Tanner graph to be a trapping set can be easily derived. 
Before describing the recursive procedure for constructing 
trapping set profiles in the next subsection, we state the 
following proposition, which gives a sufficient condition for 
the convergence of an algorithm J 7 on a Tanner graph G. 

Proposition 2: Consider decoding with an algorithm T on 
a Tanner graph G. Let V e be the set of initially corrupt variable 
nodes and / be the induced subgraph on V c . Then, algorithm 
T will converge after at most Z™ decoding iterations if there 
does not exist a subset V s of V(G) such that V s D V c and the 
induced subgraph on V s is isomorphic to a graph in 

Proof: Follows from the definition of ^/(J 7 ). ■ 
Remark: Proposition [2] only gives a sufficient condition 
because the existence of V s C V(G) which satisfies the 
above-mentioned conditions does not necessarily indicate that 
G e <§AT). 

B. Constructing a Trapping Set Profile 

The recursive procedure for constructing a trapping set pro- 
file relies on Proposition [T] which states that in the decoding 
of J 7 on a trapping set S with an inducing set /, every variable 
node is corrupt at the end of some iteration. As a result, 
one can start from the graph / and simultaneously perform 
decoding and add variable nodes to / until a trapping set is 
obtained. We now describe this procedure. 

Let assume that we are only interested in trapping sets with 
at most n max variable nodes. Consider the decoding of J 7 on 
a Tanner graph I with V(I) being the set of initially corrupt 
variable nodes. If J 7 fails on the subgraph I of I then = 
{1} and we have found the trapping set profile. If J 7 does not 
fail on the subgraph I of /, then we expand I by recursively 
adding variable nodes to / until a trapping set is found. During 
this process, we only add variable nodes that become corrupt 
at the end of a certain iteration. Consider all possible bipartite 



graphs obtained by adding one variable node to the graph I 
such that when the decoding is performed on these graphs 
with V(I) being the set of initially corrupt variable nodes, the 
newly added variable node is a corrupt variable node at the end 
of theirs? iteration, i.e., w^ v ^^ +1 £ {l w , l s }. Let Oj denote 
the set of such graphs. Take one graph in Oj and denote it 
by U. If T fails on the subgraph I of U then U G #/(T), 
otherwise we put U in a set of Tanner graphs denoted by E}. 
Repeat this procedure for other graphs in O} to obtain E}. Let 
us now consider a graph U € Ej. Again, we denote by OL the 
set of Tanner graphs obtained by adding one variable node to 
the graph U such that when the decoding is performed on these 
graphs with V(I) being the set of initially corrupt variable 
nodes, the newly added variable node is a corrupt variable node 
at the end of the first iteration, i.e., W\ v ^\ +2 £ {lw> Is}- It is 
important to note that the addition of variable node |V(I)|+2, 
which is initially correct, cannot change the fact that variable 
node | V(/)| + 1 is also corrupt at the end of the first iteration. 
This is because the addition of correct variable nodes to a 
graph does not change the states of the existing check nodes. 
We now take a graph in O- and determine if it is a trapping 
set with / being the inducing set. If it is not a trapping set 
then we say that it is a member of E|. Therefore, all graphs 
in Ej can be enumerated. In a similar fashion, we can obtain 
E|, E|, . . . , jr;^" max- ' v ^ / -''- ) . For the sake of convenience, we 
also let E° = {/}. 

At this stage, we have considered one decoding iteration on 
/. It can be seen that if S is a trapping set with at most n max 
variable nodes then either S has been found, or S must contain 
a graph in Ufej ax- '^^' ^ E}. Therefore, we proceed by 
expanding graphs in y(™™ ax_ l v '( / )l _1 ) _ 

Let K denote a Tanner graph in y("™ ax- l 1/ ( / )l~ 1 ) E}. 
We now repeat the above-mentioned process to expand K. 
Specifically, we first obtain O^, which is defined as the set 
of all Tanner graphs obtained by adding one variable node to 
the graph K such that when decoding is performed on these 
graphs with V(I) being the set of initially corrupt variable 
nodes, the newly added variable node is a corrupt variable node 
at the end of the second iteration, but not a corrupt variable 
node at the end of the first iteration. Graphs in O^- that are not 
trapping sets form the set E^ . By recursively adding variable 
nodes, graphs in ~E 2 K , E 3 K , . . . , E^" ax are enumerated. 

One can see that there are two recursive algorithms. 
The first algorithm enumerates graphs in E^- from 
graphs in EV. The second algorithm enumerates graphs in 

y(rw x -|V(4j) ^ from gj . aphs ^ y(rW-inOI-l) Jji. Each 

recursion of the second algorithm corresponds to a decoding 
algorithm. As a result, the trapping set profile is obtained after 
I™ recursions of the second algorithm. 

When expanding a graph by adding variable nodes, it is 
critical that the algorithm result in distinct graphs. In other 
words, although there can be many possible ways to add a 
new variable node, some will result in the same trapping set. 
Details will be given in the journal version of this paper. 



V. Selecting Two-Bit Bit Flipping Algorithms 

Due to page limits, we only summarize the most important 
criteria for selecting two-bit bit flipping algorithms. Let us first 
briefly discuss the number of possible algorithms. 

A. On the Number of Algorithms 

Let Q be the set of all functions from A v x 5d v — > A v that 
satisfy the symmetry and the irreducibility condition. Due the 
symmetry condition, \Q\ < 4 2x l~ d vl. There are two possible 
values of A v , and two possible values of A c . However, with 
a given A c , the two sets of algorithms T that correspond to 
two possible A v are identical (as S and W , l s and l w can 
be interchanged). Consequently, if we disregard the maximum 
number of iterations, then \&\ = 2|Q| < 2 (4|Hd vl +1 ). One can 
easily show that |S(jJ = ( dv ^ 3 ). Therefore, an upper-bound 
on the number of two-bit bit flipping algorithms is: 



For example, this upper-bound is 2 81 when d v — 3, and is 
2 141 when d v = 4. 

Due to the huge number of possible algorithms, it is 
necessary to focus on a small subset of algorithms. This subset 
of algorithms may be obtained by imposing certain constraints 
on the function /. One example of such a constraint is as 
follows: if f(0„E,) G {l w ,l s } then f(O w ,Q £ {l w ,l s }. 
This constraint requires that when a strong zero variable 
node is flipped with a given combination of check nodes, a 
weak variable node is also flipped with the same check node 
combination. Other constraints on / are derived by analyzing 
possible transitions of variable nodes and check nodes for a 
small number of iterations. 

B. Selecting a Two-Bit Bit Flipping Algorithm 

We first discuss the main criterion to select one algorithms 
among all possible algorithm. Let nf™ be the smallest number 
of variable nodes of Tanner graphs in We would like 

to select an algorithm T such that rif™ is maximized. The 
justification for this selection criterion relies on the following 
proposition, whose proof is omitted due to page limits. 

Proposition 3: Given three random Tanner graph G, Si , S2 
with < \V(Si)\ < \V(S 2 )\ < \V(G)\, the probability that 
G contains S2 is less than the probability that G contains S\. 

From Proposition [3] one can see that the larger the number 
1^(5)1 of a given Tanner graph S is, the easier it would be 
(if at all possible) to construct a Tanner graph G that does not 
contain S. Therefore, a larger rt"™ means that the sufficient 
condition for the convergence of T can be met with higher 
probability. In this sense, an algorithm T with a larger n" 1 ™ 
is more favorable. 

If for two algorithms T\ and F2, nf 1 ^ = nf^ , then 
one can derive other comparison criteria based on ^/(.Fi) 
and Sj{T2), and/or compare T\ and T2 with a different 
assumption of /. For example, the probability of a graph G 
containing a trapping set S can be also be evaluated based on 
|C(5)|. We now discuss the selection of multiple algorithms. 



C. Selecting Multiple Two-Bit Bit Flipping Algorithms 

In this paper, we only consider decoder T> with algorithms 
J~i, J-2, ■ ■ ■ , J- p operating in parallel, i.e., the received vector 
of the channel is the input vector for all algorithms. Trapping 
sets and trapping set profiles for a such decoder can be defined 
in the same manner as trapping sets and trapping set profiles 
for a two-bit bit flipping algorithm. One can show that the 
trapping set profile with an inducing set I of the decoder T> 
can be constructed with the following steps. For simplicity, 
we assume p — 2 and note that these steps can be generalized 
to include more algorithms. 

1) Construct Sj{Fi). 

2) Consider the algorithm Ti and expand each graph 
U £ S"j(Fi) using the recursive procedure described in 
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Section IV-B| with V(I) being the set of initially corrupt 
variable nodes. We then obtain a set of graphs which 
we denote by g\(T\, J-2). Algorithm Ti fails on the 
subgraph 7 of a graph in <§\{T\,T<2). 
3) Remove all graphs §\(T\, T2) of which T\ does not 
fail on the subgraph I. We then obtain the trapping set 
profile S\{T>) of the decoder T>. 
Now we can select a decoder T> among other decoders based 

on the smallest number of variable nodes of a Tanner graph 

in Sj{V). 

Remark: Knowledge on the Tanner graph of a code C can 
be used in the selection of algorithms. For example, if it is 
known that the Tanner graph of C does not contain a certain 
subgraph Y, then all graphs containing Y must be removed 
from a trapping set profile. 

VI. Numerical Results 

As an example, we describe a selection of two-bit bit 
flipping algorithms for regular column-weight-three LDPC 
codes with girth g = 8. For simplicity, we let A v = (0 S , l s ), 
A c = (0 p ,lp) and Iji = 30 for all algorithms. By imposing 
certain constraints on the functions /, we obtain a set of 
21,962,496 two-bit bit flipping algorithms. Out of these, 
there are 360, 162 algorithms which can correct any weight- 
three error pattern. Such an algorithm is capable of correcting 
any weight-three error pattern because its trapping set profile 
S'j(J') with any inducing set I containing three variable nodes 
is empty. Since all weight-three error patterns can be corrected 
with a single algorithm, our next step is to select a collection of 
algorithms which can collectively correct weight-four and -five 
error patterns with high probability. To achieve this goal, we 
construct all trapping set profiles with inducing sets containing 
four and five variable nodes for each algorithm. Note that there 
are 10 possible inducing sets (Tanner graphs with girth g — 8) 
containing four variable nodes and 24 possible inducing sets 
containing five variable nodes. Hence, for each algorithm, we 
construct a total of 34 trapping set profiles. From the trapping 
set profiles of all algorithms, we select a collection of 35 
algorithms based on the criterion mentioned in the previous 
section. Then, we simulate the performance of a decoder T> 
which employs these algorithms in parallel. The maximum 
total number of iterations of T> is 35 x 30 = 1050. 



& tn- 



10 



i 10- 6 



10 



10 




Cross-over probability (a) 

Fig. 1: Frame error rate performance of the decoder T>. 

Figure [T] shows the frame error rate (FER) performance of 
T) on the (155, 64) Tanner code. This code has d v — 3, d c = 5 
and minimum distance d m i U = 20. For comparison, the FER 
performance of the SPA with a maximum of 100 iterations 
is also included. It can be seen that the FER performance of 
V surpasses that of the SPA in the error floor region. It is 
also important to note that if we eliminate all trapping sets 
containing subgraphs that are not present in the Tanner graph 
of this code, then all the obtained trapping set profiles are 
empty. This indicates that T> can correct any error pattern up 
to weight 5 in the Tanner code. 

Figure [T] also shows the FER performance of 2? on a quasi- 
cyclic code C732 of length n = 732, rate R — 0.75 and 
minimum distance d m i n = 12. The FER performance of the 
SPA is also included for comparison. It can be seen that the 
slope of the FER curve of T> in the error floor region is higher 
than that of the SPA. This indicates that the FER performance 
of T> would eventually surpass that of the SPA. Finally, we 
remark that the slope of the FER curve of T> in the error floor 
region is between 5 and 6, which indicates that T> can correct 
error patterns of weight 4 and 5 with high probability. This 
also agrees with the fact that in our simulation, no weight-four 
error pattern that leads to decoding failure of T> was observed. 
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